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ABSTRACT 



We present the first numerical results on the binary 
coalescence of a quark star with a black hole, ob- 
tained with a 3-D Newtonian smooth particle hy- 
dro (SPH) code. The star is initially represented 
by 17,000 particles modeling a self-gravitating fluid 
with the equation of state P — [p ~ po)c^ /3, and 
the black hole by a point mass with an absorbing 
boundary at the Schwarzschild radius. As in similar 
calculations carried out for a stiff polytrope, the stel- 
lar core survives the initial episode of mass transfer, 
but here an accretion disk is clearly formed as well. 



1. INTRODUCTION 



Binary coalescence of compact objects must give rise 
to gravitational radiation emission which is expected 
to be observed with such instruments as LIGO or 
VIRGO (e.g. Thorne 1995). Neutron stars co- 
alescing with black holes may also be the source 
of r-process nuclei and of gamma-ray bursts (Lat- 
timer and Schramm 1976). Coalescence of any com- 
pact objects (black holes, neutron stars, quark stars) 
may release enough energy to power a fireball giv- 
ing rise to some cosmological gamma-ray bursts 
(Paczynski 1986), and the coalescence of strange 
(quark) stars may be especially promising (Haensel 
et al. 1991), provided that no quark matter is spilled 
into the environment with dire consequences for the 
birth rate of neutron stars (Caldwell and Fried- 
man 1991, Kluzniak 1994). 

More detailed motivation and further references can 
be found in Lee and Kluzniak (1999a), where the 
SPH coalescence of a tidally locked neutron star and 
a black hole is discussed. Lee and Kluzniak (1999b) 
discuss the coalescence of a tidally locked stiff poly- 
trope with a black hole, while SPH simulations of 
non-rotating neutron stars in binary coalescence 
with black holes have been presented by Lee (2000). 



2. NUMERICAL METHOD AND INITIAL 
CONDITIONS 



We have used the Smooth Particle Hydrodynam- 
ics (SPH) method (Monaghan 1992) for our calcula- 
tions. Our Newtonian code is three dimensional and 
we include a term in the equations of motion that 
simulates the loss of angular momentum to gravi- 
tational radiation in the binary (this is calculated 
in the quadrupole approximation for a point-mass 
binary). We also compute the gravitational radia- 
tion waveforms emitted during the coalesence, in the 
quadrupole approximation. The strange star is mod- 
eled with N = 17000 SPH particles, and has mass 
Mss = 1.4Mq and radius Rss = 11 km. The black 
hole is modeled as a point mass producing a New- 
tonian potential, and we model accretion onto it by 
placing an absorbing boundary at the Schwarzschild 
radius rsch = "^GMbh/c^- For the strange star we 
take the equation of state to be P = c^(p — po)/3, 
with po = 5 X 10^''' kg m~^. 

To perform a dynamical coalescence calculation, we 
place the binary components in a close Keplerian or- 
bit at a small separation (a few stellar radii) and give 
them the corresponding orbital velocities for point 
masses. The strange star is initially spherical and 
has no spin as seen from an external, inertial refer- 
ence frame. The orbit decays through the emission of 
gravitational waves and because of tidal instabilites 
of purely Newtonian origin (Lai, Rasio & Shapiro 
1993). 

Here we present the result of one calculation, with an 
initial mass ratio q = AIss/Mbh = 0.3 and initial 
separation = 3Rss- The dynamical evolution of 
the system is followed for approximately 7 ms. 



3. RESULTS 



At the start of the simulation a tidal bulge appears 
on the strange star as a result of the presence of 
the black hole companion. The orbital decay causes 
the star to overflow its Roche lobe, beginning mass 
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could be dynamically ejected from the system, by 

computing the total mechanical energy (kinetic + 
gravitational potential) of the fluid. We find that 
5 X 10~^Mq could leave the system, and would orig- 
inate from the outermost section of the largest blob 
formed at the end of the tidal tail. The rest of 
the matter (0.15 Mq) is gravitationally bound to 
the black hole, located practically at the center of 
mass of the system, and which has acquired a kick 
velocity from the asymmetric interaction (ffeicfe — 
lO^km s-i). 

Some of the matter in the accretion torus has enough 
angular momentum so as not to be accreted immedi- 
atly by the black hole (i.e. it can remain in orbit in 
a stable fashion), amounting to O.OSMq. This fluid 
would presumably be accreted onto the black hole 
on the viscous time scale as angular momentum is 
redistributed in the disk. 



Figure 1. Mass accretion rate onto the black hole as 
a function of time. 
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transfer within the first orbit. The star is tidally 
stretched into a tube of practically uniform thickness, 
and some matter that moves towards the black hole 
winds around it and produces a torus (see Figure 1). 
The obiter parts of the star move away from the hole 
as the stream connecting the two becomes thinner. 
Small condensations form at regular intervals along 
the narrow stream (see panel (d) in Figure 1), rem- 
iniscent of the results of a black hole neutron star 
coalescence for a stiff equation of state (Lee 2000). 

The mass accretion rate onto the black hole shows 
a large initial peak at the initial mass transfer 
episode (see Figure 2), and a subsequent decay, due 
to the presence of the accretion torus. The small 
spikes in the curve (such as the ones visible at 
t ~ 2.1,2.4,3 ms) indicate the accretion of small 
knots along the matter stream coming from the star. 

In the accretion disk around the black hole, the den- 
sity has dropped below po and so this matter behaves 
as pressure-free dust in orbit about the hole (for ref- 
erence, in Figure 1 log(/9o//0*) = —0.62). However, 
in the surviving core, flung to a large distance during 
the tidal encounter, and in the small condensations 
seen in the matter stream connecting the two binary 
components, the density remains above this value, 
and thus the matter continues to feel hydrodynamic 
effects. By t ~ 5 ms, there are 0.15 solar masses 
in the system excluding the black hole, and this is 
roughly evenly divided between matter with densi- 
ties above and below po- 
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4. DISCUSSION AND CONCLUSIONS 



At the end of the calculation we have attempted 
to obtain an estimate for the amount matter that 




Figure 2. Density contour plots in the orbital plane at (a) t=0.85 m,s; (h) t=1.27 ms; (c) t=1.7 ms; (d) 
t=2.12 ms. The axes are labeled in km. All contours are logarithmic and equally spaced every 0.5 dex. Bold 
contours are plotted at log(p/p*) = -4, -3, -2, -1 (with = Mss/Rss = 2.08 x lO^^kg m~^). 



